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Abstract 

We consider a chain described by a next-nearest-neighbor hopping combined 
with a nearest-neighbor spin flip. In two dimensions this three-body term 
arises from a mapping of the three-band Hubbard model for Cu02 planes 
to a generalized t — J model and for large 0-0 hopping favors resonance- 
valence-bond superconductivity of predominantly d-wave symmetry. Solving 
the ground state and low-energy excitations by analytical and numerical meth- 
ods we find that the chain is a Luther-Emery liquid with correlation exponent 
K p = {2- n) 2 /2, where n is the particle density. 

PACS number: 74.20.-z, 71.27.+a, 75.10. Jm 
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Recently, there has been a considerable interest in obtaining the exact ground state of 
different strongly correlated systems, using a variety of techniques Although in some 

cases the model or the parameters are rather unrealistic, exact results are useful in clarifying 
the effect of different physical ingredients and as a test of approximations. Particularly 
interesting are those models in which the ground state is superconducting, or has dominant 
superconducting correlations in one dimension @J^], because they display superconducting 
mechanisms which are different from the conventional electron-phonon one and might be 
relevant to the understanding of the high-T c systems. 

Here we consider a one-dimensional (ID) system described by the Hamiltonian: 

H t "=t" ^ 4+6>a C iSa(- - ■ Si+6), (1) 

iS^S'cr / 

where in terms of ordinary fermion operators c\ a = c| CT (l — cl_ a Ci- a ) (double occupancy at 
any site is forbidden), Sj is the spin operator at site i, and i + 8 denote the nearest neighbors 
of site i. Eq.(l) can also be written in the form of a hopping of singlet pairs: 

H t „ = 2t" J2 blb lS ; b\ 5 = -^(4 + * T 4| - cU4r) ( 2 ) 
In 2D, H t /i was obtained as an important term in a systematic analytical derivation 
of a generalized t — J model, starting from the extended three-band Hubbard model |TU 



Also, a numerical study of a CU4O8 cluster has shown that fitting of the energy levels using 
a generalized t — J model improves considerably, and a correct ordering of the levels is 



obtained, if H t /i is included [11]. If the direct 0-0 hopping t pp of the three-band model is 



large enough, t" has opposite sign as the one which corresponds to a canonical transformation 



of the Hubbard model [|1(4|T2|| and simultaneously the direct nearest-neighbor hopping t' has 
the sign which agrees with the observed Fermi surface and other properties of hole-doped 
superconductors, and with the asymmetry in the magnetic properties between hole-doped 



and electron- doped cuprates p3| , |i~2| -|l~~il . 



A recent exact diagonalization of a 4x4 cluster shows that when \t"/t\ > 0.12 with the 
sign of t" corresponding to large t pp , the system shows strong signals of <i-wave superconduc- 
tivity for realistic J/t = 0.4 and doping x = 1 — n = 0.25 [12|. The <i-wave superconducting 



susceptibility is strongly enhanced and the ground-state wave function has a significant 
overlap with a a simple short-range resonance-valence-bond state with superconducting off- 
diagonal long-range order (ODLRO) in the (i-wave and s-wave channels (the former being 
the dominant). The fact that H t " changes the physics dramatically is also evident in the 
symmetry of the ground state as a function of doping [|HJ . Another case of the three-band 
Hubbard model in wich H t " is important is the limit t pp << t p d << A where A is the 
Cu-0 charge transfer energy. In this case t" becomes of the order of t, as can be seen 
from a canonical transformation that eliminates double occupancy in the effective one-band 
Hubbard model ]15| . 



Unfortunately, although the saturation of the superconducting correlation functions in 
Ref. |HJ is an indication of ODLRO, a study of a small cluster cannot assure the presence of 



ODLRO in the thermodynamic limit. Instead, in ID one can use powerful results obtained 
using conformal-field theory to relate the correlation exponents with thermodynamic quan- 
tities and then, the effect of H t n on superconducting correlations can be determined 
unambiguosly. As a first step, we drop the other terms of the generalized t — J model. This 
allows us to obtain analytical results and to isolate the effect of H t n. 

To gain physical insight into the ID problem, let us consider a chain with only one 
vacant site (or hole). The action of the Hamiltonian Eq.(2) over a state in which the hole is 
localized at a given site, can be thought as a permutation of the hole with a nearest-neighbor 
singlet. Thus, to take maximum advantage of H t n, it is convenient that the spin configuration 
consists of a bond-ordering wave of non-overlaping singlets at each side of the hole (see Fig. 
1). The subsequent motion of the hole does not alter the spin configuration, the problem 
becomes equivalent to a tight-binding model with next-nearest-neighbor hopping 2t", and 
the resulting ground-state energy in the thermodynamic limit —At" coincides with the lower 
bound of H t n which can be obtained using Gerschgorin's theorem P7|,5||. This fact confirms 



that the spin configuration is the optimum one. When more than one hole is present in 
the system, in order to have all spins paired in singlets of nearest-neighbor particles, it is 
necessary that the distance between two successive holes is odd (the parity of this distance 



is conserved by Hf). From the arguments given above, it seems clear that the ground state 
is a combination of sequences of these singlets and holes (see Fig.l). Although we could not 
prove analytically this statement for more than two holes, we have verified that the ground- 
state energy of an open chain with 14 sites and any even number of particles is the same as 
that of a spinless model (Eqs.(3) and (4)), obtained using the mapping sketched in Fig.l: 
each singlet is replaced by a spinless fermion. For a chain of length L and N particles (N 
even) , the hopping, length and number of particles of the effective spinless fermion problem 
are given by: 

t = 2t"; L' — L — N/2; N' = N/2. (3) 

For example in Fig.l L = 14, iV = 10, and then V — 9, N' — 5. In a chain with open 
boundary conditions, the order of the singlets is conserved, and thus, the statistics of the 
fermions in the equivalent problem does not play any role. We recall that the ground-state 
energy of a spinless chain with open boundary conditions of length L', N' fermions and 
hopping t is (see for example Ref. [§): 

(2N'+1 \ 
sni7r , Jo \ , . 

I*" 1 )' (4) 

Eqs.(3) and (4) describe the ground-state energy of H t n for a chain with open boundary 
conditions. However, for the discussion below concerning the central charge and low-energy 
excitations, it is more convenient to know the energy of H t " in a ring (periodic boundary 
conditions). In this case, the mapping to the fermionic system can still be done: the fermionic 
sign which arises under a cyclic permutation of the fermions can be absorbed taking periodic 
(antiperiodic) boundary conditions if the number of fermions is odd (even). This corresponds 
to an even number of particles not multiple of four (multiple of four) in the original system 
described by H t ". In both cases the energy can be written as: 

M^/V) (5) 

1 1 sin(7r/L') v ' 

Using Eq.(3), the ground-state energy per site e = E/L of a ring of length L and density 
n = N/ L described by H t n is: 



4\t"\ sin(^) 

L S1I H L(2=n)) 

We have also verified Eq.(6) for L = 14 and all n multiple of 1/7. From Eqs. (3) and 
(4) or Eq. (6) one obtains in the thermodynamic limit: 

e(n,oo) = (2-n)sin(- ) (7) 

7i 2 — n 

The corresponding wave function can be obtained through the mapping procedure in a 
similar way as in Ref. ||. 

To obtain the long-distance behavior of any correlation function, using results of 
conformal-field theory, one needs information about the low-energy excitations of the system 
|i~6| . It is easy to see that the low-energy triplet excitations on the ground state, act as bar- 
riers to the propagation of the holes. From arguments similar to those presented above for 
the one-hole case, one sees that the spin excitation of lowest energy is to change one singlet 
of nearest-neighbor particles into a triplet. For open boundary conditions in the equivalent 
fermionic problem, this is equivalent to changing a chain of length 2L' + 1 with 2N' + 1 
particles into two equal chains of length V whith N' particles. Using Eqs. (3) and (4) we 
obtain that the value of the spin gap in the thermodynamic limit is: 

A s = |2i"|[l + cos(^^)] (8) 

The gap is maximum for n = and vanishes at n = 1. 

The charge excitations, and in particular the charge velocity v c can also be calculated us- 
ing the equivalent model of spinless fermions. It should be taken into account that distances 
and therefore also momentum are changed in the mapping procedure. While the minimum 
excitation energy E min of H t n in a ring corresponds to q = 2tt/L, the corresponding momen- 
tum for the same energy in the equivalent spinless ring is q' = 2tt/L'. Calling v' c the charge 
velocity of the spinless ring, we can write ||18|| : 



E m i n — E _ L E min — E _ L t 
2tt/L ~ U 2tt/L' ~ U Vc [ ' 
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In the thermodynamic limit v' c = 2t sm(7rN' / L'), and using Eqs.(3) and (9) one obtains: 

v c = sin 10 

2-n K 2-n J v ; 

Expanding the denominator of Eq.(6) up to third order in the argument and using Eqs.(7) 
and (10), one can see that for L — ► oo: 

e(n, L) = e(n, oo) - (11) 



where the number c = 1 is called the central charge flq]19|l . This scaling law characterizes 



a Luther- Emery liquid f20fl . In this liquid, the correlation functions are characterized by 



a single exponent K p . The superconducting pair-pair correlation functions depend on the 
distance d between pairs as d~ l ' Kp for d >> 1 and are the dominant correlation functions at 
large distances if K p > 1 [p~6| , p0| , p~8|1 . In this case a small coupling between different chains 
gives rise to superconductivity. K p can be calculated in terms of the charge velocity and the 
compressibility k flTEyigfl . Using Eq.(7) on obtains: 

1 irv c d 2 e 87i\t"\ rm 
^2 = 2K~ p = Qrf = (2-nf Sm( 2^ ) ' ( ' 

and using Eq.(10), a very simple expression for the correlation exponent results: 

K^^f- (13) 

This exact result allows to test numerical methods and approximations. As an example, 
we have calculated K p in a ring of 14 sites from the expression K p = 7r(Kn 2 Z) c /2) 1 / 2 , obtaining 



k and the Drude weight D c numerically as in Ref. |18 |. The comparison with Eq.(13) is shown 
in Fig. 2. One can see that the agreement is very good except for a small deviation at the 
extreme densities. This deviation is mainly due to the discretization of the expression of the 
second derivative which enters the compressibility. 

The system has dominant superconducting correlations for < n < 2 — \/2. Inspection 
of Eq.(12) shows that there is no phase separation in the system for any density. This is in 
contrast to the t — J and other models proposed to the cuprates in which superconductivity 
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takes place near a region of phase separation f21|j . In the dilute limit n — > 0, K p — > 2, as 
in the negative U Hubbard model for large \U\: the physics is like that of hard-core bosons 
which undergo Bose condensation. The same type of physics is expected in 2D in the dilute 
limit. We have confirmed this by numerical calculations in a 4x4 cluster and comparison 
with a BCS-like wave function. Instead, for n — > 1, K p — > 1/2, as in the case of spinless 
fermions. Addition of an attractive term like the superexchange J of the t — J model should 
extend the range of densities for which K p > 1 to larger values of n. One also expects that 
addition of J should close the spin gap for n ~ 1. The effect of the interplay of H t " with the 
terms present in the t — J model in ID has been studied quite recently by numerical methods 



| 22| , |23| . The main result is that the regions of superconductivity and phase separation for 
n ~ 1 shift to higher (lower) values of J/t if a small t" is added with the same (opposite) sign 
as the one which corresponds to a canonical transformation of the Hubbard model. Thus, 
as in 2D, the sign of t" favored for large t pp in the mapping from the three-band Hubbard 



model fllQ| , |12|| is more convenient for superconductivity for realistic doping. 
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FIGURE CAPTIONS 



Fig. l:(a) Schematic representation of a typical state with definite occupation at each site 
among those which compose the ground state of the system, (b) state obtained from (a) 
replacing each singlet pair by a spinless fermion. Open circles denote vacant sites, closed 
circles joined by a full line represent a singlet b\ s (Eq. (2)), and up arrows represent spinless 
fermions. 

Fig. 2: Correlation exponent K p obtained numerically in a ring of 14 sites (solid squares) in 
comparison with the exact result (solid line). 
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(b) o f f oo t t t o 



Fig. 1 
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Fig. 2 



